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We study interacting dipolar atomic bosons in a triple- well potential within a ring geometry. This 
system is shown to be equivalent to a three-site Bose-Hubbard model. We analyze the ground state 
of dipolar bosons by varying the effective on-site interaction. This analysis is performed both nu- 
merically and analytically by using suitable coherent-state representations of the ground state. The 
latter exhibits a variety of forms ranging from the su(3) coherent state in the delocalization regime 
to a macroscopic cat-like state with fully localized populations, passing for a coexistence regime 
where the ground state displays a mixed character. We characterize the quantum correlations 
of the ground state from the bi-partition perspective. We calculate both numerically and analyti- 
cally (within the previous coherent-state representation) the single-site entanglement entropy which, 
among various interesting properties, exhibits a maximum value in correspondence to the transition 
from the cat-like to the coexistence regime. In the latter case, we show that the ground-state mixed 
form corresponds, semiclassically, to an energy exhibiting two almost-degenerate minima. 

PACS numbers: 03.75.Ss,03.75.Hh,64.75.+g 



I. INTRODUCTION 

The physics of quantum gases is dramatically affected 
by intcratom.ic interactions [l|. These, in most cases, are 
effectively described by a short-range isotropic potential 
characterized by the s-wave scattering length. However, 
when the quantum gases are dipolar, then the anisotropic 
dipole-dipole interaction between magnetic or electric 
dipole moments becomes very important due to its long- 
range character and determines novel interesting proper- 
ties i-i. 

Very recently, dipolar bosons trapped by triple- well po- 
tentials have became a topic of growing interest. This 
system often features the open-chain geometry where all 
condensates are mutually coupled through dipole inter- 
action but only the two lateral condensates exchange 
bosons with the central one. Among many interesting 
aspects, the triple-well system exhibits states revealing 
the non-local character of the the dipole-dipole interac- 
tion within the Josephson-like dynamics |5| and a com- 
plex ground-state phase diagram depending on inter-site 
and on-site interaction |6| . 

A different interesting choice for the space structure 
of the triple well is the ring geometry characterizing the 
closed chain with periodic boundary conditions. Its re- 
alization by means of optical lattices has been designed 
in Ref. [3] and is supported by the recent construction 
of traps able to confine bosons in toroidal domains g . 
Such a geometry has been widely used for investigating 
the so-called Bose-Hubbard (BH) triple-well where only 
on-site (nondipolar) interaction are present. The lattice 
geometry reflects on the dynamics which exhibits a rich 
variety of translational invariant states and includes ring 
currents represented by vortex states [9J-[14|. The ring 



geometry, however, gives its most dramatic manifestation 
in the attractive strong-interaction regime where, owing 
to the inherent translation symmetry, the ground-state 
is a Schrodinger cat formed by a superposition of macro- 
scopic localized states [ij, [T^] . Recently, the transition 
between different density-wave phases has been studied 
for dipolar bosons in a ring-shaped lattice [13J . 

The ring geometry of the model representing dipolar 
bosons plays a crucial role in establishing its complete 
equivalence with the symmetric BH triple-well. Simple 
calculations show how, with M = 3 sites, the term repre- 
senting dipolar interaction can be absorbed by the on-site 
interaction term. This process can be effected neither 
with M > 3 nor with an open chain. This equivalence 
gives the possibility to exploit the considerable amount of 
information about the dynamics of symmetric triple- well 
to investigate three-well dipolar models. 

In this work, we consider a number N of interacting 
dipolar bosons at zero temperature confined by a triple- 
well potential forming an equilateral triangle. The mi- 
croscopic dynamics of this system is effectively described 
by a 3-site extended Bose-Hubbard (EBH) Hamiltonian 
which includes the nearest-neighbor hopping J and the 
density-density interaction on the same site t/o and be- 
tween different sites C/i as well. We diagonalize this 
Hamiltonian and study the ground-state properties by 
varying the effective on-site interaction U = Uq — Ui. 
When this is zero, the ground state of the EBH Hamilto- 
nian reduces to a su(3) coherent state (describing the bo- 
son delocalization) which tends to a Fock state (with the 
same boson number in each well) when U > becomes 
large enough. We show that, for positive U, the ground 
state exhibits the well-known Mott-like (superfluid) form 
when the effective interaction is sufficiently large (weak). 



For negative U, the ground state goes from a su(3) coher- 
ent state (where bosons are delocahzed) to a macroscopic 
superposition of three locahzed states (cat-like state) 
which, for \U\ large enough, evolves towards the NNN 
state (superposition of three fully-localized Fock states, 
representing the generalization of the well-known NOON 
state). In the intermediate U range between such states, 
the ground state is expected to have a mixed character 
[14| . We systematically explore this range determining its 
extension and making fully visible the complex form of 
the ground-state in which a coherent state (describing a 
uniform boson distribution) is superimposed to a macro- 
scopic cat-like state. The semiclassical analysis of the 
low-energy scenario in the cat-like regime reveals that the 
latter features the coexistence of two minima which can 
be identified as the classical counterpart of the ground 
state with a mixed character. 

The variety of ground states previously described can 
be characterized from a genuine quantum-correlation 
point of view. We shall calculate the single-site entan- 
glement entropy S, both numerically, by diagonalizing 
exactly the EBH, and analytically (in the three ground- 
state regimes: the NNN state, the coexistence, and the 
su(3) coherent state) in the limit of large N b y fo llowing 
the same approach of QQ, [l3|- In Refs. J17l - [l9| it was 
shown that the Schrodinger cat state (NOON) is not the 
most bi-partite entangled state among the BH Hamil- 
tonian ground-states. Likewise, we point out that the 
NNN-like state does not exhibit the maximum bi-partite 
entanglement entropy: S, actually, attains its maximum 
value within the coexistence regime. We have analyzed as 
well the first derivative of S with respect to the interac- 
tion by relating its minimum and maximum, respectively, 
to the emergence of three Fock localized states (around 
the delocahzed state) and the precursor of NNN state. 



tact potential T4r(r — r') ~ gS{Y — r') and a long-range 
dipole-dipole potential 



Vm{v - r') = 7 



1 - 3 cos^ e 
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(2) 



Here g = Ai:h^as/m with a^ the interatomic s-wave scat- 
tering length, 7 = fiQ^^/An for magnetic dipoles (^o is 
the vacuum magnetic susceptibility and /x is the mag- 
netic dipole moment) or 7 = d'^/Aneo for electric dipoles 
{eq is the vacuum dielectric constant and d is the electric 
dipole moment). For sufficiently large external (electric 
or magnetic) fields the boson dipoles are aligned along 
the same direction, so that 9 is the angle between the 
vector r — r' and the dipole orientation. 

The model describing dipolar bosons can be derived 
from the bosonic-field Hamiltonian 

H = I d^r¥{r)Ho^{r) 

dhd^r'^^r)^''{r')V(r - r')*(r')*(r) 

Ho^^—V^+Vtrapir), (3) 

where operators ^(r) and ^'''(r) annihilates and creates 
a boson at point r, respectively. Operator ^(r) can be 
expanded in terms of space modes 

^(r) = 01 (r) ai + (j)2 (r) as + 03 (r) 03 , (4) 

where space modes dk, a|. satisfy the Heisenberg-Weyl 
algebra [dk, d'U = Skq while, due to the form of the trap- 
ping potential, single-particle wave functions (/)fc(r) take 
the form 



II. THE MODEL HAMILTONIAN 

We consider a system of A^ dipolar dilute interacting 
bosons trapped by a potential Vtrap(r) confining bosons 
in a three- well array. Vtrap (r) is obtained by superposing 
a strong harmonic confinement along axis z with three 
(planar) potential wells located at the vertices of an equi- 
lateral triangle with side £ in the x ~ y plane. Each of 
these wells is described by the product of two Gaussians 
with the same width w. Thus Vtrap{'^) reads 



ytrapV^) 



,2.^2 



Vb ^ cxp 



2{r-nf 



(1) 



with m the mass of each boson, uJz the trapping frequency 
in the axial direction, Vo the depth of each well and ri — 
(£/2,0), ra = (0,73^/2) and r^ = (-^/2,0). We focus 
on the case w <^ £ entailing strong localization in the 
proximity of sites r^ . 

The boson-boson interaction potential V{r — r') = 
Vsrir — r') -I- Vd(j(r — r') is the sum of a short-range con- 



0fe(i") = 9{z)wk{x,y) 



(5) 



in which g[z) represents the ground-state wave func- 
tion of harmonic potential {rwjj1/2)z'^, and Wk{x,y) is 
a single-particle wave function localized in the kth well. 
Functions Wk{x,y) satisfy the orthonormality condition 
j dxdyw^{x,y)wi{x,y) = Ski since they are suitable lin- 
ear combinations of the three lowest energy eigenfunc- 
tions of the (planar part of) triple- well potential Vtrap- 
Hence / d3r(/)*(r)0,(r) = Ski- 

By assuming symmetric wells, the resulting dipolar- 
boson model is described by the effective 3-site extended 
Bose-Hubbard (EBH) Hamiltonian 



H = —J[d\d2 



Ssfti 



aids 



0302 



dlds + agSi] 



U, 



"9" X! "'(""^ - 1) + C/i [nih2 + h2m + him] ■ (6) 



In equation ^ fik — 0^,0^ counts the number of particles 
in the fcth well. 



Microscopic processes are described in Hamiltonian (|6]) 
by the three macroscopic parameters J, Uq and Ui. The 
hopping amphtude J is given by 



J 



dh^kir) 



|^V2 + Vtrapir) 



Mr) , (7) 



where k and I [k ^l) arc two nearest-neighbor sites. 

Concerning the on-site interaction C/q, the latter em- 
bodies the contributions of both short-range and dipole- 
dipole interactions (notice that in a single cloud the 
dipole-dipole interaction yields a mean-field energy shift 
which is strongly dependent on the anisotropy of this in- 
teraction [5|)- Then, [/q is given by 

+ 7 /d3rdV|0fc(r)pFdd(r-r')|0fc(r')P ■ (8) 



The third parameter is the nearest-neighbor interaction 
Ui. In principle there would be an exponentially small 
contribution of the short-range potential since the over- 
lap between the wave functions on two adjacent wells is 
nonzero. This, however, can be shown to be completely 
negligible @ . Then one finds 



Ui^^Jd^rdh'\Mr)\''Vdd{r- 



r')\Mr) 



(9) 



Interestingly, owing to the underlying ring geometry, 
the dipolar-boson model 1^ can be mapped onto a sim- 
pler model in which dipolar-interaction terms nifik are 
completely removed. Since 



iV^ 



3 



rij + 2{hih2 + n2n3 + hifi^) 



where iV = fii +fi2 -(-na commutes with Hamiltonian ([B]), 
then H reduces to a BH model 

H = — J[a{a2 + dgfli + a^^'i + <^3'^2 + d\aj, + d\ai\ 



^J2n.,{n,-l)+^N{N-l), (10) 



where U = Uo — Ui shows that the nearest-neighbor 
(dipolar) interactions have been absorbed by the on-site 
interaction term. In the following, we will neglect the 
constant term UiN{N — l)/2. Then, the boson-boson 
interaction can be described in terms of an effective on- 
site interaction of amplitude Ua ~ Ui. The latter is re- 
pulsive when the Uq > Ui and attractive in the oppo- 
site case. This means (see the discussion below) that 
making Uo larger or smaller than Ui allows one to ex- 
plore a variety of ground-states ranging from the Mott- 
like state (for U ^ J /N and integer filling) and attrac- 
tive/repulsive superfluid states {\U\ < J/N), involving 
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FIG. 1: (Color online). Coefficients Icna.nsI s-s functions 
of occupation numbers 712 and 713 for different values of r = 
J/{N\U\) and A'^ = 45 bosons. Top-bottom (eacfi row, from 
left to rigfit). First row (repulsive case U > 0): r 3> 1, 
r = 0.000222. Second row: r = 0.2415, r = 0.2389. Third 
row: r = 0.2374, r = 0.2349. Fourth row: r = 0.2315, r = 
0.2020. The last three rows concern the attractive case U < 0. 



the same average population per site, to Schrodinger-cat 
states (|C/| ^ J/N with U < 0) consisting essentially 
of the superposition of three fully-localized macroscopic 
states. 



III. DIPOLAR-BOSON GROUND STATE 

The ground-state of the Hamiltonian (jTU]) is a super- 
position of different Fock states which, due to the con- 
servation of total boson number N, can be written as 

I*) = E E ^"2."3l^-("2+'^3),'^2,n3) . (11) 
n2— ns— 

State 1^) features different regimes that depend on inter- 
action parameters Uq and Ui through the effective hop- 
ping amphtute r = J/{N\U\). 



For T ^ 1, namely ii \U\ <IC J/N (note that Uq can ap- 
proach Ui either from above or from below), the ground 
state tends to the su(3) coherent state [2^ 



|del) 



1 



/M 



T. 



N 



z,a]] |0,0,0), (12) 



where Zi — l/v3, and |0,0, 0) is the vacuum state with 
zero bosons. State |del) describes the delocalization of 
boson population since the boson-number expectation 
value {z\ni\z) = iVjZip = A^/3 is the same at each well. 
Note that in the limit \U\ — > 0, state p^ represents an 
exact eigenstate of PH]) . In this case, in fact, Hamilto- 
nian H can be rewritten a,s H = —J{3A^A — N) with 

A — {a[ + 0,2 + a\)/y3 which, owing to its harmonic- 
oscillator form, is diagonalized by the state (|12p . 

When r ^ 1 and U < the ground state of Hamilto- 
nian (fTUl) tends to a linear combination of the Fock states 
each one describing the complete localization of bosons 
at a single well. Then |^) ~ \W), where 

\W) - i=(|7V,o,0) + |0,iV,0) + |0,0,7V)) . (13) 

State \W), hereafter called NNN state, can be interpreted 
as the A^-boson generalization of the three-qubit state 
[13 \W) = {\1,0,0) + \0,1,0) + \0,0,1))/V^. Since each 
component of \W) features a macroscopic boson localiza- 
tion, such a state represents a Schrodinger cat. 

Finally, in the presence of integer filling, for repulsive 
interaction [/ > and r <C 1, ground state (fTTj) takes the 
form 



1^) ~ |n, n,n) , n 



N 



eNo 



(14) 



where the majority component \n,n,n) reflects the bo- 
son space distribution characterizing the Mott phase. For 
T = one has |^) = \n,n,n). If the filling is non integer 
the ground state is well represented by the symmetric 
form 1^) ~ {\n + s,n,n) + \7i,n + s,n) + \n,n, n + s))/\/3. 
This describes the expected superfluid character through 
the presence of a fully-delocalized particle (hole) excita- 
tion when s = +1 {s = —1). 

Fig. [T] shows the squared modulus |c„2_„3p of the coef- 
ficients of 1^) components as a function of occupation 
numbers n2, and 713. The total number of bosons is 
TV = 45 entailing integer filling. The left top panel of 
Fig. [TJ where U = describes an infinitely large r, rep- 
resents the su(3) coherent state (|12|) . Starting from this 
configuration, the smaller r the narrower the distribu- 



tion of \c„ 



In general, for a generic r S [0,oo] in 



the repulsive regime C/ > 0, the quantity |c„2_„3p attains 
its maximum value for n2 = n^ = 15. This situation is 
evidenced in the right top panel for r = 0.00022, where 
the Mott-like ground state is very close to state ([T4|. 

For attractive interactions U < and r ^ 1 the 
ground state is, in general, a delocalized state qualita- 
tively similar to that of case U > with a finite r 3> 1 



(see the left top panel of Fig. [T]). Below the critical value 
To (one has tq = 0.25 for N large enough) ground state 
(111!) totally modifies its structure exhibiting a transition 
[14| . |15| to a mixed state characterized by the coexis- 
tence of three symmetric localized states and a delocal- 
ized state (well represented by state (HI])). 

When T is further decreased (top-bottom, second- 
fourth rows of Fig. [1]), ground state PT|) once more 
changes its structure tending to a Schrodinger-cat state 
essentially coinciding with state p^ . Translated into 
the language of |c„2_„3p, the left panel of the second row 
(r = 0.2415) well illustrates the emergence of three lo- 
calized states: the distribution of |c„2,„3p, in fact, starts 
to show three lateral peaks in addition to the central 
one describing a uniform boson distribution. By further 
decreasing interaction [/ < (sec Fig. [1] second row, 
r = 0.2389) the emergence of three external peaks be- 
comes more and more evident. 

The third row of Fig. [1] clearly shows how the stronger 
the interatomic attraction \U\, the smaller the central 
peak. In our sequence, the latter reaches its smallest 
nonzero value for r = 0.2349 (third row, right panel) and 
then disappears for r = 0.2315 (fourth row, left panel) 
entailing a ground state characterized by three well sepa- 
rated peaks. When r becomes small enough, the ground 
state evolves towards a form well represented by the NNN 
state (see the right bottom panel where t = 0.2020). 

Some considerations are now in order concerning effec- 
tive interaction U = Uq — Ui. First, it is worth noting 
that for any value of U the same ground state (and more 
in general any eigenstate of H) can be realized by (in 
principle) infinitely-many choices of Uq and Ui if such 
parameters are changed by the same amount. Interest- 
ingly, while in the BH model (where Ui = 0) the rich 
structure of the attractive ground state is obtained if in- 
teraction Uq is genuinely negative, in the EBH model, 
due to dipolar interactions, the condition U < can be 
met with both Uq and Ui positive. A nontrivial conse- 
quence is that one can investigate the attractive regimes 
of triple-well dynamics avoiding the possible collapse of 
the bosonic cloud expected, for sufficiently large atomic 
densities, in the presence of truly attractive bosons. 



IV. ENTANGLEMENT PROPERTIES OF I*) 

The dramatic changes characterizing the structure of 
the dipolar-boson ground state when r is varied suggest 
that quantum correlations between different sites may 
play an important role. To this end we analyze the entan- 
glement properties pertaining to state ([TT|) for different 
choices of interaction parameters. The density matrix p 
associated to 1^) is 



/5=l*)(*l 



(15) 



We describe the quantum correlations by calculating the 
entanglement between the well 1 and the system formed 



by the wells 2 and 3. An excellent measure of this en- 
tanglement is provided by the single-site entanglement 
(SSE) entropy S [22|- This quantity is the von Neumann 
entropy of reduced density matrix pi defined by 



pi = Tr2Tr^p 



(16) 



The latter is a matrix obtained by partially tracing the 
total density matrix ([T5|) . first over the degrees of freedom 
of the well 3 and then on those of the well 2. Therefore, 
the SSE entropy S is given by 



5" = -Tri(/5ilog2/5i) 



(17) 



We have studied S" as a function of t both for attrac- 
tive (Fig. [21 panel a) and for repulsive (Fig. [21 panel b) 
interactions. For repulsive bosons (Fig. [21 panel b), S 
increases when r is increased. S tends to an asymptotic 
value (not shown) confirming that the ground state tends, 
in parallel, for r ^ 1 to state |del) with a nonzero, con- 
stant quantum-correlation degree reflecting the nonsepa- 
rablc, Gaussian form of \^) in the superfluid regime. For 
T — )■ the vanishing of S reflects the disentangled form 
of Mott-like state \n,n^n) described by a single integer- 
filling Fock state. 

For attractive bosons, S exhibits three remarkably 
points showed in Fig. [21 panel c. For t = 0.2416, the 
(first) r derivative of S has a minimum. This value of r 
is very close to that for which the three external peaks 
crop up (Fig. [H second row, left panel) around the cen- 
tral one inherited from the regime with r 3> 1. When 
T = 0.2360, the single-site entanglement entropy attains 
its maximum value {S = 5.0201). In correspondence to 
this value the central peak of the distribution of |c„2^„3 p 
reaches its smallest nonzero value (Fig. [1] third row, right 
panel) . Finally, the first derivative of S has a maximum 
for T = 0.2306. This value is very close to the value of 
T for which the NNN state precursor emerges, that is 
the state where the central peak disappears (see Fig. [I] 
fourth row, left panel). Then, the NNN state is achieved 
by further lowering r. 

It is possible to provide analytical estimations of the 
SSE entropy in the limit of large N . To this end we 
consider the generic su(3) coherent state 



l^) = 



1 



/M 



(Ciol + ^20; 



■iAr\ 



(18) 



where Sf^^l^ip = 1 ensures the normalization condition 
(^I'l^I') = 1. Within the present coherent-state picture 
the expectation value Ni = (\l/|ni|^) — N\^i\'^ oi fii can 
be interpreted as the ith-site population. By defining the 
relative phases 0ij = 6i — 6j 



J'Pij 



6g 



(19) 



the expectation value E = ('^\H\'^) of Hamiltonian ^TU\\ 
on the state (TlSt is found to be [201 



U 



E^-N{N-1)J2\^, 






■04*) (20) 
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X 

FIG. 2: Top panels. The entanglement entropy S (|17p vs. 
T — J/{N\U\) for N=4:5 bosons. Panel a): attractive interac- 
tion U < 0. Panel b): repulsive interaction L*' > 0. Panels c): 
the first derivative of S vs r in the attractive case L'^ < 0. 



whose final expression in terms of variables Nk and (jtij 

is the well-known boson Josephson-j unction form E = 

u 

2 



1 



■ E. N! - 2JZk<, yAWV^cos(<^,,.) if {N - 1)/N 
is assumed. Based on ansatz (fT8)) wc evaluate the SSE en- 
tropy ([T7| of the ground state for three significant cases: 
the (uniform) ground state described by state (fT2l) . the 
cat-like ground state (including the case of the NNN state 
([rsl)). the regime with the coexistence of different ground 
states. 



V. ANALYTIC ESTIMATION OF ENTROPY S 
IN DIFFERENT GROUND-STATE REGIMES 



To begin with, we derive the analytic expression of the 
reduced density matrix p£ when |^) is a coherent state 
(|T8| . This expression is obtained by partitioning the sites 
in two blocks [1, £] and [.^-1-1,3], with £ = 1, 2 and tracing 
out the sites belonging to the second block [Ifl, [I4I 



- E 



n- 



nf+i,...ii3 1.1.1=1+1 ''^i- 



\D{n,+^,...n^)W (21) 



with the constraint X]i=i "-i = -^ ^-nd where 



D(n,+i,...n3) = (0| n (aO"'l*)- 

The resulting pf^ is & m x m diagonal matrix with m = 
{N + ^)!/(^!7V!). In particular, for £ = 1, the diagonal 
elements of this matrix are found to be jl6| 

Pi,«(ei) = (^)(i-ian''^""^iar (22) 

This expression gives the probability to find n bosons 
in the well 1, and satisfies the normalization condition 

Sn=o/°i." = 1- I'^ ^^^ large- A'' limit this binomial dis- 
tribution can be approximated to a Gaussian one with 
mean value -/V|CiP and variance -/V|^ip|(l — |^ip). For 
N ^f CO, therefore, the entropy can be calculated analyt- 
ically. In this limit, it exhibits the following asymptotic 
expression [l6| 



5^1log2(27reiV|aP(l-|6P) 



A. Delocalized-boson ground state 



(23) 



For C/ = (r — > oo), the ground state is found by 
setting $1=^2=^3 = l/%/3 in formula (fT8|) . This gives 
simply 



|del) = 



1 



^t , -t , -UN 



v¥nr\ 



{a\+al + air\0), (24) 



which is state ([T^ describing the boson delocalization. 
From formula (P^ . with |^ip = 1/3 and large iV, we find 
that the entanglement entropy in the delocalized-boson 
regime is given by 



Sdei ^ 2 l°g2 (47reiV) - loga 3. 



(25) 



This value exactly reproduces the asymptotic value ob- 
tained numerically for S when t ^ 1. 



B. Cat-like ground state 

The strong- interaction regime (|(7| ^ J 1^) is charac- 
terized by a ground state representing a cat-like state. 
The latter, tending to the NNN state dTS]) for r ^ 0, can 
be constructed by exploiting state (|18p . Classically, for 
T small enough, the self-trapping effect entails an almost 
complete localization of attractive bosons at some site 
(soliton-like state). Quantum-mechanically, the peak of 
the boson population can be placed, with the same prob- 
ability, in one of the three sites of ring. A realistic ansatz 
for this state thus consists in superposing three states 



([TS]) each one involving the localization at a different site. 
One has 



I cat) 



1 



vsivT 



E(^4+^(4+i+4+2))'^|0), (26) 



fc=i 



where 0^+3 = a^ within the periodic ring geometry. The 

fcth state |*fc) = (^4 + ^("Ui+4-H2))^|0)/^^ofthis 
superposition features local populations 



(*fc|nfe|*i 



m? 



(*i.|n,;|*fc) 



nH- 



for i = fc + l,fc -f- 2. In each |^fc), well k exhibits a 
population larger than the those of wells A: -I- 1 and k + 2, 
namely \iu? = \i? > l^+iP = 1^-^21' = \v\^ where Cj 
are the coherent-state parameters of definition p^ . 

For TV ^ 1 the coherent-state product reduces to the 
simple form {'^hl'^k) — ^hk while in the energy expecta- 
tion value 

£;= (cat|H|cat) = i^^(vI/,|F|*,,) 



mixed terms [^ k\H\^ h) with h ^ k can be shown [1J| 
to vanish for large N. In addition, due to the translation 
simmetry oi H, {'^k\H\'^k) yields the same expression 
for each k. Then, for fc = 1, one finds E = (cat|i^|cat) 
~ (*i|iJ|*i) giving 



E: 



^N{N 



imf + 2\rf)~2JN{^r^*+Cri+\rf)- 



Interestingly, the latter coincides with the energy ob- 
tained from formula (|20p when setting £,1 = £, and 
^3 = ^2 = V- This makes it evident that the effective en- 
ergy E is essentially independent from the symmetrized 



form of cat-like state ((26 
and the inverse formulas 



By defining z ~ \£Y 



m 



\e 



2z + l 



(27) 



one obtains 
U 



Ec^-N{N~l) 
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l + 2z' 



A^(l-z + 2v/(l-z)(l + 2z)cos(/)i2) , (28) 



2J 
~"3~ 

where cf>i2 = (f>i — (f>2 is defined by formula (|19p . Note that 
only such a relative phase survives in E since ^2 = Cs — ^ 
in l^i). In order to identify the energy minimum one 
must impose 012 = 0, phase ^12 being independent from 
z. The stationary values of energy ([28]) in the attractive 
case U < 0, for iV ^ 1. are found by the solutions of 



f->^ 



t(42 - 1) 



ynr^Zy/l + 2z 



= 0, (29) 



with r = J/(\U\N). For r < 1, in addition to z = 
0, equation ([^^ exhibits two solutions one of which is 



negative. In this regime, these are necessarily placed in 
the proximity of asymptotes z = 1 and z = —1/2 of 
function dE/dz. The negative solution zi ~ —1/2 can be 
discarded in that, owing to |^p ^ jyyp, it corresponds to 
an antisoliton state (this typically represents an excited 
state) where one of the three wells is essentially depleted. 
In this regime solution, z = (which satisfies equation 
(j29|l for any r) and zi correspond to an energy maximum 
and to a local minimum, respectively. This can be easily 
checked through the calculation of the second derivative 
E = d^E/dz^. The positive solution can be found by 
substituting z = I — v (with u <C 1) in equation (|29p . 
This gives 



Z = 22 ~ 1 



3r2 



(i-ry 



\e^l-2T\ \v\'^r\ 



which can be shown to represent the effective minimum 
of E{z) being E > 0. The fact that the energy minimum 
entails |.Jp ^ jr^p is consistent with the distribution of 
bosons prescribed by cat-like state ((26)) . This validates 
the ground-state representation based on formula (|26| . 

The entanglement entropy is found by assuming that 
the reduced density matrix, in the large N limit, is 
formed by a symmetric linear combination of three 
almost-normal distributions, i.e. 



/5cat = :r(pi(a)+/5i(6) + /5i(6)) 



(30) 



(matrix pi is defined by formula ([H])) refiecting the pa- 
rameter choice ^1 = £,, ^2 = ^3 ^ V of coherent state I*!*!). 
In the extreme case of the NNN state, namely for z = 1, 
the entropy is equal to (logj 3 — 2/3). Interestingly, the 
same result is achieved by calculating the single-site en- 
tropy 5 = — J2n=a Pi,nlog2Pi,n ^^ ^hc State NNN where 
pifi ~ 2/3, pi^N = 1/3 and pi^„ = otherwise. This 
seems to confirm the validity of the form assumed for 
the density matrix pcat- More in general, for a generic 
cat-state in the large N limit, we can write 



Scat ^ log2 3 - 2/3 + 5(0/3 + 25(77)/3 , 



(31) 



where S{a) — — Tri (/5i(a) log2 /5i(a)), with a ~ ^ 77, 
are the single-site entropies for sites populated by iV|^p 
and A^|?7p bosons, respectively. In the limit z — > 1 these 
entropies are zero since they corresponds to Fock states. 
Hence we recover the result for the NNN state. Instead, 
as soon as 1 — z 3> 1/A'^ (where z, however, must ensure 
a population |^p considerably larger than |?7p). S{£^) and 
S{ii) are given by Eq. (P5)) with £,1 = £, and ^ = 77. In 
this case Eq. (|3T|) becomes 



Jl0g2 
Jl0g2 



47reA^(l-z)(l + 2z) 
27reA^(l-z)(2 + z) 



(32) 



C. Semiclassical interpretation of the transition to 
the cat-like state 



Equation ((29| . derived from energy pS]). supplies fur- 
ther interesting information about the transition from 
the cat-like to the uniform state for t -^ 1/4. If r is 
increased, the negative solution zi approaches solution 
z — and for r = 2/9 the two solutions coincide. This 
value is particularly significant because solutions zi and 
z = change their character when zi crosses z = 0. 

For r > 2/9 solution zi becomes positive while Z2 < 1 
moves away from 1. By assuming that zi is positive but 
small, equation (|29p supplies the approximate expression 
zi = 3(t - 2/9)/r giving £;(zi) < 0. Then zi > cor- 
responds to an energy maximum. In parallel, a second 
minimum crops up at z = 0, in addition to the energy 
minimum associated to solution Z2. At z = the second 
derivative is found to be ^(0) = -2\U\N^{1 - 9r/2)/3 
confirming that this solution becomes a minimum if 
T > 2/9. The relevant populations, given by formulas 
(|27p with z = 0, show how the second minimum involves 
the uniform distribution corresponding to state (|12p. 

The occurrence of a second energy minimum when pa- 
rameter T is varied is extremely interesting. It is the sig- 
nal that, within the present coherent-state picture of |^), 
the system manifests a second way to realize a minimum- 
energy configuration in the range 2/9 < r < 1/4. In this 
interval the numerical values of E{0) and E[z2) are very 
close. For r = 1/4 the exact solutions of equation (^51) 
can be easily found. These are zi = 1/4 and Z2 = 1/2 in 
addition to z = 0. In this special case £^(0) = E{z2). 

This behavior, emerging in the interval [2/9, 1/4], can 
be viewed as the semiclassical counterpart of the consid- 
erable changes characterizing the ground-state structure 
shown in Fig. [1] for 0.2020 < t < 0.2415. The pres- 
ence of two (semiclassical) energy minima corresponds, 
quantum-mechanically, to the coexistence in the ground 
state of two dominating components noted in Ref . [ij] : 
the cat-like state (with three peaks describing the three 
possible ways to get a complete localizations) and state 
P^ characterized by the boson delocalization. 



D. Coexistence of |cat) and |del) 

The previous discussion in subsection IV CI and the re- 
sults illustrated in Fig. [1] suggest that, in the regime clas- 
sically identified by 2/9 < t < 1/4, the ground-state fea- 
tures the coexistence of a cat state and a (dclocalized) 
coherent state. Explicitly, our ansatz for the ground state 
I*) is 



|mix) = V 1 — Q^^ |cat) + a |del) 



(33) 



where state |del) is given by ((T2|) and a = l/-\/2 has 
been assumed. Moreover, we consider the state |cat) cor- 
responding to solution z — Z2 = 1/2 of case r = 1/4, 
where parameters £ and 7/ are such that |^p = 2/3 and 
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|77p = 1/6. By construction the overlap between state 
I cat) and the locahzed coherent state |del) is small and 
become negligible in the limit of large N. This follows 
when the formula describing the product of two coher- 
ent states (^IC) = {J2i=i^iCi) is applied to the product 
of I del) with one of the three components of |cat), and 
among the three components of |cat) themselves. Under 
this assumption, the entanglement entropy is given by 

Smtx ^ (1 - a^) Scatiz = 1/2) + a^Sdel 

-(1 - a2) log2(l - a2) - a^ loga a^ , (34) 



where Scat{z = 1/2) is given by Eq. ([5^ by setting 
z = 1/2 and Sdei is given by Eq. ^. For a^ = 1/2, Eq. 
([M| becomes simply 

5™„ ~ i log2 (47reiV/3) + ^ (l + logs 5) • (35) 

Notice that Smix is always greater than the entropy of 
the cat state, Eg. ([5^ . for any value of z. For N = 45, 
Smix — 5.05, which is in a quantitatively good agreement 
with the maximum value achieved by 5, found through 
the exact diagonalization of the BH Hamiltonian (see the 
comments about Fig. [2]). 



VI. CONCLUSIONS 

We have studied a system of interacting dipolar bosons 
confined by a triple-well potential with periodic bound- 
ary conditions. The latter can be described in terms of 
a 3-site extended Bose-Hubbard (BH) model including 
the boson interaction Ui between nearest-neighbor sites. 
Thanks to its symmetry properties, we have reduced the 
dipolar-boson model to the symmetric BH model depend- 
ing on a unique (effective) interaction U. Experimentally, 
this fact is certainly interesting because one can explore 
the attractive regime {U < 0) of triple- well dynamics 
by using positive-valued onsite and nearest-neighbor in- 
teractions (C/q and t/i, respectively). This circumstance 
allows one to study the attractive regime of BH model 
avoiding the possible collapse due to Uq < for suffi- 
ciently large atomic densities. 

We have analyzed the ground-state structure of 
dipolar-bosons in a triple-well by changing the (unique) 



parameter U in its whole variation range. Fig. [T] illus- 
trates the numerical study of the ground-state structure 
for different U. In the range [/ > 0, the ground state 
has been shown to exhibit the expected Mott-insulator 
(superfluid) form for r = J/(JV|C/|) < 1 (r > 1). For at- 
tractive interaction [/ < 0, the ground state goes from a 
su(3) coherent state (describing the boson delocalization 
of the superfluid state) to a macroscopic superposition of 
three states with localized populations (the cat-like state) 
separated by an intermediate range where the ground 
state has a mixed character. The latter is well described 
by a linear combination of the delocalized-boson ground 
state and the macroscopic cat-like state. The semiclas- 
sical study of the low-energy landscape in the cat-like 
regime shows that the energy features the coexistence of 
two minima which can be viewed as the classical coun- 
terpart of the ground state with a mixed character. 

To stress the importance of the quantum correlations 
between different sites, we have investigated the single- 
site entanglement entropy as a function of the effective 
interaction. This analysis has been carried out both nu- 
merically, by diagonalizing the extended BH Hamilto- 
nian, and analytically, by representing the ground state 
in terms of coherent states. In particular, we have found 
that the single-site entanglement entropy S reaches its 
maximum within the coexistence regime. The S maxi- 
mum is attained in correspondence to the transition from 
the four-peak to the three-peak structure of the ground 
state. In parallel, the r derivative S is found to have a 
maximum and a minimum point in the vicinity of which 
S undergoes a rapid change. The rapid increasing of S 
takes place in the regime heralding the NNN state, while 
the rapid decreasing of S corresponds to the merging of 
the lateral peaks with the central one. 
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